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This  report  s ummarize s the  research  carried  out  during  the  past  three 
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l 4-7  5-C-0392  | The  details  of  the  experimental 
work  can  be  found  in  the  Wayne  State  University  Ph.D.  Dissertation  of  Wolody- 


myr  Czubatyj , 1977  , Department  of  Electrical  and  Computer  Engineering.  The 
important  conclusions,  both  experimental  and  theoretical,  are  contained  in  the 
attached  reprints,  which  comprise  the  major  portion  of  this  report. 

In  brief,  we  experimentally  investigated  transferred  electron  effects  in 
bulk  n-GaAs  and  n-InP  under  hydrostatic  pressures,  and  compared  the  results 
with  the  results  of  our  numerical  Monte  Carlo  calculations.  Our  experimental 
results  differed  substantially  from  those  of  other  investigators,  and  our  numeri- 
cal results  showed  that  the  discrepancies  were  due  to  the  presence  of  substan- 
tial ionized  impurity  scattering  in  our  samples.  In  general,  we  not  only  verified 
that  fundamental  band  structure  effects  are  masked  by  contact  phenomena,  but 
also  that  they  are  seriously  affected  by  ionized  impurity  scattering. 

To  obtain  useful  results  for  fundamental  band  structure  parameters  using 
transferred  electron  effects,  the  investigator  must  begin  with  high  purity,  high 
mobility  material,  relatively  free  from  ionized  impurities.  This  can  usually  be 
achieved  by  epitaxial  growth.  The  epitaxial  samples  should  be  grown  on  a semi- 
insulating  substrate  and  shaped  into  *H*  or  turret  patterns  to  minimize  effects  . . 
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"due  to  the  metallic  contacts.  Turret  shaped  samples  are  best  since  in  one 
polarity  a boundary  condition  producing  either  transit  time  oscillations  or 
current  saturation  due  to  a high  cathode  field  will  yield  values  of  the  satur- 
ated drift  velocity.  In  the  opposite  polarity,  the  absence  of  a domain  mode 
and  a sufficiently  high  value  of  the  peak  velocity,  along  with  current  satur- 
ation or  switching,  or  circuit  controlled  oscillations,  can  then  provide  con- 
fidence that  low  boundary  field  conditions  have  been  achieved. 
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TRANSFERRED  ELECTRON  EFFECTS  IN  n-GaAs  AND 
/i-InP  UNDER  HYDROSTATIC  PRESSUREt 


W.  Czubatyj.  M.  S.  Shur  and  M.  P.  Shaw 

Department  of  Electrical  and  Computer  Engineering.  Wayne  Stale  University,  Detroit.  Ml  48202.  U S A. 

Abstract  -Experimental  results  are  presented  showing  the  variation  of  the  peak  velocity,  electric  held  at  peak 
velocity  and  saturated  velocity  of  the  velocity-held  curve  for  hulk  n-GaAs  under  hydrostatic  pressure.  Similar  data 
are  presented  for  the  held  at  peak  velocity  for  hulk  n-lnP.  The  extraction  of  the  data  from  the  current-voltage 
characteristics  is  based  on  an  exploitation  of  the  influence  of  boundary  conditions  on  the  manifestation  of 
transferred-electron  induced  (Gunn)  current  instabilities.  Our  results,  which  differ  from  those  of  other  in- 
vestigations. are  compared  with  our  Monte  Carlo  calculations  of  the  velocity-held  curves  for  these  materials.  We 
hnd  that  ionized  impurity  scattering  plays  an  important  role  in  the  understanding  of  our  data 


1.  INTRODUCTION 

Two  recent  experimental  investigations  of  the  high  elec- 
tric field  transport  properties  of  ri-GaAs| 1. 2]  and  n- 
InP[3]  under  hydrostatic  pressure  have  helped  provide 
further  clarification  of  the  position  and  variation  with 
pressure  of  some  conduction  band  extrema  in  these 
materials.  Furthermore.  Vinson  el  a/.[l]  have  also  made 
uniaxial  stress  measurements  and  obtained  numerical 
evidence  supporting  the  view  that  the  lowest  subsidiary 
minimum  in  the  conduction  band  of  GaAs  lies  at  the  L 
point  in  the  Brillouin  zone[4).  In  the  above  investigation, 
the  high  quality  material  employed  was  epitaxially 
grown  from  the  vapor  phase  onto  it'  substrates  and  the 
top  contact  was  an  alloyed  metal  (Sn  or  Ag-Sn)  uni- 
formly covering  the  entire  active  region.  In  our  studies 
we  have  used  less  pure,  long  bulk  (horizontal  Bridgman 
and  Czochralski)  samples  where  contact  effects  have 
been  minimized  by  geometrically  shaping  the  sample  to 
remove  the  active  region  from  the  influence  of  (he 
cathode  contact.  Although  both  sets  of  results  exhibit 
some  of  the  same  qualitative  features,  there  are  sufficient 
substantial  quantitative  differences  so  that  it  is  difficult 
to  obtain  an  unambiguous  interpretation  of  the  data. 
These  differences  will  be  emphasized  in  the  following 
two  sections  where  our  experimental  results  are  presen- 
ted In  the  last  section  we  discuss  the  results  of  our 
numerical  calculations,  which  demonstrate  that  ionized 
impurity  scattering  plays  an  important  role  in  controlling 
the  character  of  our  experimental  results. 

2.  EXPERIMENTAL  RESULTS  FOR 
The  experimental  studies  involve  applied  voltages 
sufficiently  high  to  induce  the  Gunn  instability (5).  Vinson 
et  al. |l|  and  Pickering  el  a/.|2]  have  both  observed  that 
the  threshold  voltage.  V„  for  the  Gunn  instability  in- 
creases with  applied  hydrostatic  pressure,  P,  in  the  range 
from  0 up  to  about  15  kbar.  At  higher  values  of  P,  V,  first 
decreases  then  increases  with  P until  the  instability 
disappears  near  30  kbar.  The  observation  of  a range 
where  V,  decreases  with  P is  consistent  with  the  original 
observations  of  Hutson  et  a/.|6].  It  was  demonstrated  by 
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Pickering  el  al.  that  the  failure  to  observe  an  initial 
increase  in  V,  with  P by  Hutson  et  al.  was  due  lo  (he 
nonlinear  type  of  alloyed  metal-semiconductor  used. 
Recently,  we  pointed  out  |7]  that  even  with  the  linear 
metal-semiconductor  contacts  employed  by  Vinson  et  al. 
and  Pickering  et  al..  major  features  of  the  velocity-field, 
HE).  characteristics  may  still  be  masked.  We  found  that 
once  the  controlling  influence  of  the  contacts  were  eli- 
minated. the  field  at  peak  velocity,  Er,  (see  Fig.  I)  was 
found  to  increase  substantially  with  P up  to  near 
10  kbar.  Furthermore,  the  peak  velocity,  vp.  and 
saturated  velocity.  i\,  also  increased  as  well. 

Our  experimental  and  numerical  results,  as  well  as 
those  of  Vinson  et  al.,  are  summarized  in  Fig.  I.  To 
obtain  our  experimental  results,  it  was  necessary  to 
shape  (sculpture)  the  sample  as  shown  in  the  inset.  The 
reason  for  this  is  that  unless  the  cathode  boundary 
contact  fields  are  reduced  below  Ep[ 5,7],  a transii-time 
instability  will  occur  at  currents  below  the  current  den- 
sity jr  associated  with  vp.  Samples  that  were  not  sculp- 
tured produced  either  current  saturation  or  transit-time 
effects. 

Our  specimens  were  cut  from  previously  metallized 
and  alloyed  wafers  of  | 111]  oriented  bulk  material  having 
carrier  concentrations  in  the  range  1.7-2. 3 x 10"  cm  ' 
and  mobilities  in  the  range  3000-5500  cm2/V-scc.  All  our 
samples  exhibited  an  exponential  increase  in  resistance 
with  P.  beginning  above  about  I kbar.  which  we  attri- 
buted to  carrier  freeze-out  to  impurity  levels  as  the  band 
gap  increased.  This  phenomenon  makes  the  data  analysis 
more  difficult  and  the  results  somewhat  less  precise,  but, 
as  we  shall  show,  the  results  were  consistent  with  our 
assumptions. 

The  contacts  were  formed  by  first  evaporating  onto  the 
samples  a thin  coal  of  Sn(  - 500A  thick),  then  clec- 
trolessly  depositing  Ni(  ~ I000A  thick),  and  finally 
evaporating  another  layer  of  Sn(  ~ I2.000A  thick).  The 
samples  were  then  alloyed  for  typically  three  minutes  at 
500°C.  They  were  diced  into  parallelepipeds  approximately 
0.10x0. 10x0.04  cm  on  edge.  Finally,  each  parallelepiped 
was  sculptured  into  the  form  shown  in  the  inset  of  Fig.  I by 
using  a wire  string  saw  coated  with  abrasive  grit.  When  the 
smaller  contact  was  the  cathode  the  current  saturated  at  a 
value  j,  = nev,  |5),  where  ne  is  the  charge  density.  When 
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Pig  I.  Normalized  peak  velocity  (dots)  and  saturated  velocity  (straight  lines)  vs  normalized  electric  field  al  peak 
velocity  as  a function  of  P for  GaAs.  The  data  were  extracted  from  /l  V ) curves  for  three  samples  taken  with 
100  nsec  voltage  pulses  al  a prf  of  50  hz.  For  a typical  sample  r,  = I 97  * 10’ cm/sec  and  h'r  = 4800V/cm.  The 
crosses  beginning  at  the  P - 0 point  represent  an  average  of  the  data  appearing  in  Ref  I The  curved  solid  line 
represents  an  average  of  the  numerical  results  we  obtained  and  those  of  Ref  1;  N,  = 0 The  dashed  line  shows  our 
numerical  results  for  N,  = It)17  cm  ' The  dot-dashed  line  shows  our  numerical  results  for  = It)"1  cin  ' The  inset 
displays  the  shape  of  a typical  sample,  drawn  to  scale.  The  pre-sculptured  dimensions  are  given  in  the  text. 

Sculpturing  typically  reduced  the  cross-sectional  area  by  about  a factor  of  7. 


the  larger  contact  was  the  cathode  (here  the  active  region  is 
effectively  removed  from  the  influence  of  the  cathode 
contact),  the  /(V)  characteristics  were  linear  up  to  vr. 

After  the  /( VO  measurements  were  made,  the  samples 
were  mounted  by  their  leads  on  a transistor  header, 
immersed  in  epoxy,  capped  and  centrifuged  to  remove 
any  bubbles  present  in  the  epoxy.  They  were  then  placed 
in  a core  drilled  from  the  center  of  a cut  pyrophylite 
regular  tetrahedron.  2.54  cm  on  edge.  The  core  voids 
were  filled  with  a Duco-rouge  mixture  and  the  pyramid 
rouged  and  placed  into  a tetrahedral  press  capable  of 
achieving  pressures  up  to  40  kbar.  The  system  was  cali- 
brated and  checked  for  hydrostatic  behavior  by  replacing 
the  GaAs  samples  with  Bi,  Ce.  Hg.  or  Yb.  and  recording 
the  P -induced  phase  transitions  in  these  materials. 

The  /(  V)  curves  for  both  polarities  were  measured  at 
I kbar  intervals  up  to  1 2 kbar  and  the  data  analysis 
proceeded  as  follows.  First,  at  atmospheric  P,  the  low- 
voltage  resistance.  R.  of  the  sample.  V,.  Ip.  and  /,  were 
measured,  where  1„  and  /,  are  the  peak  and  saturated 
currents,  respectively  The  ratio  /,,//,  yielded  vjc, 
directly.  vp  was  obtained  by  using  the  value  v,  = 
0.86  x 10’ cm/sec.  This  fixes  E,,(5.7|  and  hence  the 
mobility  pf=ipIEr)  The  carrier  density,  «.  was  then 
determined  from  the  resistivity,  p,  measured  prior  to 
sculpturing.  was  calculated  ( = nevp)  and  the  effective 
cross-sectional  area  A of  the  active  region  was  then 
determined  l-lpljp).  The  pscudo-samplc-lcngth  / was 
taken  as  V,IEP.  (Both  / and  A arc  expected  to  remain 
constant  to  within  1-2%  up  to  12  kbar.) 

Next,  as  P was  increased  we  determined  pi  = RAIl), 
n (using  the  fact  that  p will  decrease  by  about  10% 
from  0 to  10  kbar. [8|  which  we  checked  experimentally 
for  self-consistency),  jr.  t>„  = j,Jne).  Ep. 

Ep(  = vjp)  and  p'f  = v'pIEp).  The  primes  denote  an 


alternate  method  of  obtaining  these  parameters.  Com- 
parison of  rp  with  ip.  Ep  with  E'p  and  p with  p'  provided 
a critical  check  of  our  technique.  The  agreement  was 
never  worse  than  15%  and  typically  good  to  within  5%. 
Each  datum  point  plotted  in  Fig.  I is  an  average  of  the 
primed  and  unprimed  data.  The  sources  of  error  are 
uncertainties  in  changes  of  physical  dimensions  with 
pressure  (2%).  P(5%),  / and  V'  (3%).  and  p (10%).  It  is 
evident  that  vp.  Er.  and  r,  all  increase  substantially  with 
P up  to  10  kbar.  It  is  evident  that  ip.  Ep.  and  r,  all 
increase  substantially  with  P up  to  10  kbar.  In  the  data 
shown  in  Fig.  1 the  turnover  point  (where  |ep.  £p)  starts 
decreasing  with  P)  is  near  10  kbar.  However,  the  high  R 
of  the  sample  at  these  P's  make  the  turnover  point 
difficult  to  determine  accurately. 

J.  EXPERIMENTAL  RESULTS  FOR  n-lnP 
Pitt  and  Vyas(3|  found  that  V,  for  the  Gunn  instability 
in  InP  remained  constant  for  pressures  up  to  near 
40  kbar  After  that  V',  increased  with  P up  to  about 
70  kbar.  Our  results,  which  arc  startlingly  different,  arc 
shown  in  Fig.  2.  To  obtain  these  results,  the  InP  samples 
were  prepared  almost  exactly  as  the  GaAs  samples. 
However,  the  ''turret”  type  samples  shown  in  the  inset 
of  Fig.  I were  extremely  hard  to  fabricate  since  InP  is 
very  fragile  and  cracks  under  the  slightest  stress  along  the 
(100)  plane.  Even  when  a sample  was  successfully  fabri- 
cated. it  was  found  that  it  could  not  be  cycled  electrically 
because  the  high  fields  often  present  at  the  small  contact 
caused  irreversible  damage  Therefore,  a different  design 
was  tried  in  order  to  reduce  the  fields  al  bolh  electrodes. 
Two  deep  cuts  were  made  into  the  sample  with  one 
surface  contract  removed  (see  inset  of  Fig.  2).  (Such  a 
large  contact  area-double  cut  sample  proved  to  be  stable 
under  hydrostatic  pressure.  Voltage  probing  of  Tcledel- 
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Kg.  2 Normalized  peak  velocity  vs  normalized  electric  field  at 
peak  velocity  as  a function  of  P for  InP  The  inset  displays  the 
shape  of  a typical  sample,  drawn  to  scale. 

los  resistance  paper  cut  to  model  the  cross-sectional 
shape  of  an  actual  sample  showed  that  approximately 
half  of  the  applied  voltage  was  dropped  in  the  central 
region  and  the  remainder  in  the  bulk.)  When  the  applied 
bias  was  increased,  the  voltage  drop  in  the  central  region 
would  correspondingly  increase.  Due  to  the  large  area  of 
the  contact  region,  only  a small  voltage  was  dropped 
there.  Only  at  very  large  voltages,  approximately  two  or 
three  times  larger  than  the  peak  voltages  used,  does  the 
cathode  drop  become  significant  (25  to  33%  of  the  sam- 
ple voltage)  Upon  reaching  the  peak  current  for  the 
sample,  the  central  region  began  to  avalanche  rapidly, 
producing  in  a short  time  a sample  of  lower  resistance, 
and  yielding  an  increased  current  and  smaller  voltage 
drop  across  the  central  region.  If  the  fields  at  the  con- 
tacts were  relatively  small  at  the  time  of  this  switching 
effect,  the  switching  occurred  along  the  load  line  to  a 
stable,  higher  current  point  If.  however,  the  contact 
fields  were  too  large  when  the  switch  occurred  and  the 
redistributed  voltage  caused  breakdown  in  these  regions, 
irreversible  damage  would  occur.  Several  samples  were 
destroyed  in  this  manner 

The  data  taken  for  the  sculptured  InP  samples  (Fig.  2) 
were  similar  to  those  of  GaAs,  except  that  the  "current 
saturation  polarity"  was  absent  because  the  samples 
were  symmetric.  Kven  in  unsculptured  samples  that 
could  have  shown  saturation  for  lower  carrier  concen- 
trations, impact  ionization  within  the  high  cathode  field 
region  always  caused  the  current  to  rise  continuously 
with  increasing  voltage,  making  determination  of  a 
saturated  current  impossible  The  curves  also  had  a 
"soft"  knee,  making  determination  of  the  knee  point  an 
unreliable  process  In  general,  unsculpturcd  hulk  InP 
samples  (n22*IO,,cm  ’)  develop  sufficiently  high 
cathode  fields  so  that  the  current  observed  on  a sampling 
oscilloscope  will  usually  increase  as  the  voltage  in- 
creases. even  when  a domain  instability  is  generated. 

Our  samples  had  carrier  concentrations  typically  in  the 
range  2-4*  10"  cm  ' and  mobilities  between  3500  and 


4500  citT/V -sec.  The  resistivity  varied  from  about  0.4  to 
O.ft  11-em.  In  order  to  avoid  cumulative  breakdown 
effects  in  these  samples,  we  applied  single-shut  voltage 
pulses.  As  the  applied  bias  brought  the  current  to  near 
the  peak  value,  the  It  VO  curves  developed  a small 
sublinearity,  which  is  expected  for  InP  since  its  e(£) 
curve  is  sublincar  prior  to  threshold  (GaAs  is  linear).  At 
the  peak  current,  a switch  to  higher  current  occurs  along 
the  load  line.  We  attribute  this  effect  to  impact  ionization 
within  the  anode  high  field  region |5|. 

The  peak  current  density  J,,  = l,JA.  where  A is  known 
to  within  10%  and  /,,  to  within  3%.  yielded  a value  for  n. 
assuming  a given  value  for  tr.  This  also  yielded  a cor- 
responding p from  the  r(£)  curve  for  InP  obtained  by 
Fawcett  and  Herbert  |9|.  The  resistivity  was  then 
measured,  and  the  process  iterated  until  good  agreement 
between  r,„  p n and  p was  obtained. 

As  for  GaAs.  p as  a function  of  P was  obtained.  Next. 
JJne  yielded  cr  as  a function  of  P.  Er  as  a function  of  P 
was  obtained  by  dividing  V',  by  a constant  determined  by 
the  width  of  the  cut  and  a geometric  sculpturing  factor. 
The  data  presented  in  Fig.  2 represents  an  average  of 
three  samples. 

4.  NUMERICAL  CALCULATIONS  AND  DISCUSSION 

In  order  to  explain  the  above  data,  we  performed 
Monte  Carlo  calculations  of  the  v (£)  characteristics.  In 
this  paper  we  will  emphasize  our  results  for  GaAs.  We 
used  essentially  the  same  parameters  used  in  Ref.  I: 
Effective  mass-ni  * = 0.68m, „ mf  = 0.35m„.  m*  = 
0.35 m„;  Energy  separation-Af,  , = 0.40eV.  A Er  , = 
0.38  eV:  Pressure  dependence  of  the  energy  separations, 
(d£r_„/dP)  = - 1 1 x 10"'" eV  bar-1,  (d£r_,/dP)  = 
-5 x |()  "eV  bar  Deformation  potential  -7eV; 
Coupling  constants— the  unscreened  values  of  Ref.  9; 
Pressure  dependence  of  effective  mass.  (dnti*/dP)  from 
l p theory  with  tdEJdP)  = lOx  10  '"cVbar  1 and  m*. 
m,  constant.  (£„  is  the  energy  gap  and  I",  x , /.  are  major 
symmetry  points  in  the  Brilloiun  zone.)  Using  these 
parameters  with  parabolic  bands,  and  considering  acous- 
tic. optical,  intcrvalley  and  intravallcy  scattering,  we 
were  able  to  reproduce  the  numerical  results  of  Ref.  I 
reasonably  well  (see  Fig.  I).  Only  when  we  included 
ionized  impurity  scattering  at  a density  greater  than 
about  N,  = I0l7cm  1 were  we  able  to  obtain  results  that 
were  similar  to  ours;  an  initial  increase  in  t,,  and  a 
relatively  large  initial  increase  in  £,,  with  increasing  P 
(see  Fig.  I).  (The  use  of  the  band  parameters  suggested 
by  Aspncs|4|  made  the  agreement  between  experiment 
and  numerical  calculation  more  disparate.)  Note  that  civ,, 
remains  above  unity  over  a larger  range  of  P for  larger 
values  of  N,.  EIEr  is  also  greater  at  a given  P for  larger 
values  of  N,. 

We  suggest  the  following  model  in  order  to  explain  the 
role  that  ionized  impurity  scattering  plays  in  controlling 
the  variation  of  the  (tr,  £,,)  point  with  P.  First,  we  know 
that  En  increases  as  p decreases.  Furthermore,  since  p 
decreases  with  increasing  P,  (d£r/dP)  will  be  greater  for 
smaller  p in  the  range  where  Ep  increases  with  P.  Since 
p decreases  as  N,  increases,  we  expect  that  both  £,,  and 
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its  variation  with  P to  increase  with  increasing  N,. 

Next,  in  order  to  explain  why  ir  increases  with  in- 
creasing P.  we  assume  that  the  average  drift  velocity  is 
given  by  r = (n,  vL  + n„vm)Hn,  +«„),  where  v, . n,  and 
r„,  n„  are  the  drift  velocities  and  electron  concentrations 
in  the  lower  and  upper  valleys,  respectively.  Since  nega- 
tive differential  mobility  is  a “tail  of  the  distribution 
effect ",  then  we  expect  that  near  threshold  nu  *<  n, : also 
»„  <4  t|  Hence,  near  threshold,  v =*  v,  n,J(n,  + n„). 

We  know  that  impurity  scattering  decreases  v, . Fur- 
thermore. the  value  of  rtu  near  threshold  will  be  rather 
insensitive  to  impurity  scattering,  which  is  most  im- 
portant at  low  electron  energies [10],  Thus,  at  atmos- 
pheric P we  expect  and  observe  a lower  value  of  the 
peak  velocity  when  impurity  scattering  is  present  than  in 
its  absence  [5,  ll|.  At  low  P our  calculations  show  that  P 
increases  Ep.  with  or  without  ionized  impurity  scattering. 
At  high  fields  (near  threshold),  the  difference  between  v, 
in  the  presence  of.  and  without  ionized  impurity  scatter- 
ing becomes  smaller  (in  the  limit  of  very  high  fields  both 
values  should  be  equal).  This  will  result  in  an  initial 
increase  of  vr  with  P in  the  presence  of  impurity  scat- 
tering. The  maximum  value  of  vp  reached  under  applied 
P will  of  course  be  below  that  corresponding  to  the 
situation  without  impurity  scattering. 

In  summary,  we  have  verified  one,  and  uncovered 
another  phenomenon  which  tends  to  mask  fundamental 
band  structure  effects  in  negative  differential  mobility 
(NDM)  semiconductors.  The  new  effect  is  that  of  ionized 
impurity  scattering.  Ionized  impurity  scattering  affects 
the  peak  velocity  and  field  at  peak  velocity,  making 
accurate  determinations  of  band  structure  variations  with 
pressure  difficult  to  achieve.  The  second  effect  is  that  of 
the  influence  of  the  contacts  on  the  form  and  stability  of 
the  current-voltage  characteristics  of  NDM  semiconduc- 
tors subjected  to  hydrostatic  pressure  |2. 3, 7). 

To  improve  on  the  results  obtained  for  GaAs  and  InP, 
we  should  start  with  much  purer,  higher  mobility 


material,  relatively  free  of  ionized  impurities.  This  can  he 
achieved  by  epitaxial  growth.  The  expitaxial  samples 
should  be  grown  on  a scmi-insulaling  substrate  and 
shaped  into  "H"  or  turret  patterns  to  minimize  (or  eli- 
minate) the  metallic  contact  effects.  The  turret  shaped 
samples  (inset.  Fig.  I)  are  most  useful  since  in  one 
polarity  a boundary  condition  producing  either  transit- 
time  oscillations  or  current  saturation  due  to  a high 
cathode  field  will  yield  values  of  v,.  In  the  opposite 
polarity,  the  absence  of  a domain  mode  and  a sufficiently 
high  value  of  r„,  along  with  current  saturation  or  switch- 
ing. nr  circuit  controlled  oscillations,  can  then  provide 
confidence  that  low  boundary  field  conditions  have 
actually  been  achieved))] 
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Experimental  results  show  that  the  peak  velocity,  peak  electric  field,  and  saturated  velocity  of  the 
velocity-electric-field  curve  for  n-GaAs  all  increase  substantially  with  increasing  hydrostatic  pressure  up 
to  about  10  kbar,  the  peak  velocity  and  field  then  begin  decreasing  with  further  increase  in  pressure  The 
extraction  of  the  data  from  the  current-voltage  characteristics  is  based  on  an  exploitation  of  the  influence 
of  the  boundary  conditions  on  the  manifestation  of  transferred-electron-induced  (Gunn)  current 
instabilities. 

PACS  numbers:  72.20.Ht,  72.80.Ey 


Vinson  cl  nl.1  have  presented  numerically  supported 
experimental  evidence,  based  on  hydrostatic  pressure 
and  uniaxial  stress  measurements  at  applied  voltages 
sufficiently  high  to  induce  the  Gunn  instability,2  in  har- 
mony with  the  recent  finding  that  the  lowest  subsidiary 
minimum  in  the  conduction  band  of  GaAs  lies  at  the  L 
rather  than  the  Y point  in  the  Brillouin  zone.3  In  this 
paper  we  present  further  experimental  evidence  relating 
to  this  view  and  also  provide  additional  data  which  for 
the  first  time  displays  major  features  of  the  variation 
of  the  carrier  velocity  v versus  electric  field  E char- 
acteristic at  hydrostatic  pressures  P below  12  kbar. 
Vinson  cl  nl . and  Pickering  cl  nl  .*  have  both  observed 
that  the  threshold  voltage  VT  for  the  Gunn  instability 
increases  with  increasing  hydrostatic  P in  the  range 
from  0 up  to  about  15  kbar.  At  higher  values  of  P , VT 
first  decreases  then  increases  with  P until  the  instabil- 
ity disappears  near  30  kbar.  The  observation  of  a range 


FIG.  1.  The  vlE)  curve  and  computer  generated  j-vs-E  (aver- 
age field)  curves  (A,  D)t  IF,,  C)  for  various  values  of  Ec.  Ecl A) 

0.  £C(D,)  r,  kV/cm,  F.c (H?)  8. 5 kV/cm,  and  Ke(C)  24  kV/cm. 
Curve  A results  for  low  Ec,  curve  C for  high  F.c,  and  curves 
If,  and  ltj  for  intermediate  Ec  (in  the  NDM  region;  shaded. 

When  Ec  is  in  the  shaded  region,  the  cathodc-to-, anode  transit- 
time  instability  occurs).  The  reason  curve  A does  not  extend 
up  to  v,  is  that  small  statistical  doping  fluctuations  were  In- 
cluded in  the  calculations.  The  left-  and  right-hand  ordinates 
are  related  by  ) n^evUs),  where  n„  is  the  doping  density  and  e 
the  electronic  charge.  This  figure  is  taken  from  Refs.  2 and  0, 
where  it  was  demonstrated  that  the  calculations,  based  on  a 
fixed  Ec  model,  provide  excellent  agreement  with  experimental 
results  on  long  samples  of  n-GaAs. 
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where  VT  decreases  with  P is  consistent  with  the  orig- 
inal observations  of  Hutson  cl  nl  ,5  It  was  demonstrated 
by  Pickering  et  nl . that  the  failure  to  observe  an  initial 
increase  in  VT  with  P by  Hutson  cl  nl . was  due  to  the 
nonlinear  type  of  alloyed  metal -semiconductor  contacts 
used.  A major  point  of  our  paper  is  that  even  with  the 
linear  alloyed  metal -semiconductor  contacts  employed 
by  Vinson  cl  nl . and  Pickering  cl  nl  the  major  features 
of  the  r(E)  characteristics  are  still  masked.  Our  re- 
sults show  that  once  the  influence  of  the  contacts’5  is  (i) 
clbtiittnlcrl  and  (ii)  controlled , not  only  is  the  threshold 
field  Ep  observed  to  increase  substantially  with  P up  to 
near  10  kbar  but,  surprisingly,  the  threshold  velocity 
cp  and  saturated  velocity  r,  (see  Fig.  1)  also  increase 
as  well. 

Crucial  to  our  arguments  is  the  fact  that  the  classic 
Gunn  instability,  i.e. , the  cathode-to-anode  transit- 
time mode,3  does  not  occur  at  a critical  value  of  voltage 
nr  electric  field.  In  fact,  trying  to  extract  useful  infor- 
mation from  a variation  in  VT  for  this  type  of  instabil- 
ity is  a futile  exercise.  The  instability  occurs  at  a 
critical  value  of  current  density  jc, 2,6,1  which  is  con- 
trolled by  the  value  of  the  electric  field  in  the  negative 
differential  mobility  (NDM)  region  at  the  cathode  bound- 
ary Fc.  As  shown  in  Fig.  1,  if  Fc  is  too  high  (above  the 
NDM  region)  the  transit-time  instability  is  suppressed 
and  the  current  saturates  at  the  value  = nact\\  a high- 
field  region  remains  at  the  cathode.  In  order  to  reach 
a critical  current  equal  to  jp  (from  which  rp  can  be  ob- 
tained), low-field  boundary  conditions  must  be  achieved. 
However,  the  cathode-to-anode  transit-time  instability 
is  also  al  tost  always  suppressed  for  low  F.c\  here  a 
high-field  region  often  evolves  at  the  anode  and  the  cur- 
rent saturates  very  near  ip  = nvct'p  (Fig.  1). 

We  have  exploited  the  latter  two  boundary  condi- 
tions2'* (conlrollint;  and  clitninalinK,  respectively,  the 
influence  of  the  contacts)  in  order  to  obtain  vp,  Fp , and 
v,  as  functions  of  P.  For  our  purposes  the  most  impor- 
tant of  these  boundary  conditions  is  that  of  low  Ec.  We 
have  only  rarely  found  an  alloyed  metal -semiconductor 
contact  (for  carrier  concentrations  in  the  mid-1015  cm'5 
range  and  below)  that  produced  what  could  be  charac- 
terized as  a low  F.c  contact.  We  were  only  able  to  reli- 
ably achieve  this  boundary  condition  by  sculpturing  the 
sample  to  remove  the  active  region  from  the  influence 
of  the  contacts.*  Samples  that  were  not  sculptured  pro- 
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duced  either  current  saturation  or  transit -time  effects; 
transit -time-type  samples  (intermediate  Fc,  within  the 
NDM  region)  become  electrically  unstable  before  rp 
is  approached  and  hence  cannot  yield  information  con- 
cerning the  peak  point  (rp,  Fp).  This  is  an  important 
feature.  Guetin  and  Schreder"  showed  that  the  height  of 
a Schottky  barrier  on  n-GaAs  increases  with  increasing 
P,  implying  that  Ec  increases.  Thus,  we  expect  F.c  to 
increase  with  P for  intermediate  Fc  samples.  A larger 
Ec  implies  a lower  jc,  which  means  that  the  critical 
current  density  for  such  samples  will  decrease  with 
increasing  P while  1',-can  increase.  Such  behavior  has 
been  emphasized  in  Refs.  1 and  3,  and  plays  a major 
role  in  their  analyses. 

Our  specimens  were  cut  from  previously  metallized 
and  alloyed  wafers  of  [lllj-oriented  bulk  material  hav- 
ing carrier  concentrations  in  the  range  3000  — 5500  cm2/ 
Vsec.  All  our  samples  exhibited  an  exponential  in- 
crease in  resistance  with  P,  beginning  above  about  1 
kbar,  which  we  attributed  to  carrier  freeze-out  to  im- 
purity levels  as  the  band  gap  increased.  This  phenome- 
non makes  the  data  analysis  more  difficult  and  the  re- 
sults somewhat  less  precise,  but,  as  we  shall  show, 
the  results  were  consistent  with  our  assumptions. 

The  contacts  were  formed  by  first  evaporating  onto 
the  samples  a thin  coat  of  Sn  (-  500  A thick),  then  elec- 
trolessly  depositing  Ni  (-1000  A thick),  and  finally 
evaporating  another  layer  of  Sn  (~  12  000  A thick).  The 
samples  were  then  alloyed  for  typically  180  sec  at 
500  C.  They  were  diced  into  parallelopipeds  approxi- 
mately 0.  1 « 0. 1 x0.  04  cm  on  edge.  Finally,  each 
parallelopiped  was  sculptured  into  the  form  shown  in 
the  inset  of  Fig.  2 by  using  a wire  string  saw  coated 
with  abrasive  grit.  When  the  smaller  contact  was  the 
cathode,  the  current -voltage  l(V)  characteristics  were 
like  those  of  Fig.  1,  curve  C (high  Ec).  When  the  larger 
contact  was  the  cathode  (here  the  active  region  is  ef- 
fectively removed  from  the  influence  of  the  cathode  con- 
tact), the  l(V)  characteristics  were  those  of  Fig.  1, 
curve  A (low  Ec). 

After  the  l(V)  measurements  were  made,  the  samples 
were  mounted  by  their  leads  on  a transistor  header, 
immersed  in  epoxy,  capped,  and  centrifuged  to  remove 
any  bubbles  present  in  the  epoxy.  They  were  then 
placed  in  a core  drilled  from  the  center  of  a cut  pyro- 
phylite  regular  tetrahedron,  2.54  cm  on  edge.  The  core 
voids  were  filled  with  a Duco-rouge  mixture  and  the 
pyramid  rouged  and  placed  into  a tetrahedral  press 
capable  of  achieving  pressures  up  to  40  kbar.  The  sys- 
tem was  calibrated  and  checked  for  hydrostatic  behavior 
by  replacing  the  GaAs  samples  with  Bi,  Ce,  Hg,  or  Yb, 
and  recording  the  P-induced  phase  transitions  in  these 
materials. 

The  /(V)  curves  for  both  polarities  were  measured  at 
1 kbar  intervals  up  to  12  kbar  and  the  data  analysis 
proceeded  as  follows.  First,  at  atmospheric  P,  the 
low-voltage  resistance  R of  the  sample,  lp,  VT,  and  /, 
were  measured.  The  ratio  Ijl,  yielded  vp/v,  directly. 
rt  was  obtained  by  using  the  value  r,  = 0.  86  x 101  cm/ 
sec.  From  the  Butcher- Fawcett  r(F.)  curve  (Fig.  1) 
this  fixes  Ep  and  hence  the  mobility  p (=r^/Ef).  The 
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carrier  density  n was  then  determined  from  the  resis- 
tivity p measured  prior  to  sculpturing.  ip  was  calculated 
(=ncrp)  and  the  effective  cross-sectional  area  A of  the 
active  region  was  then  determined  (=lp/ip).  The  pseudo- 
sample-length I was  taken  as  VT/FP.  (Both  I and  A are 
expected  to  remain  constant  to  within  1—  2f|  up  to  12 
kbar.) 

Next,  as  P was  increased  we  determined  p (=  RA/I), 
it  (using  the  fact  that  p will  decrease  by  about  10'J?  from 
0 to  10  kbar,9  which  we  checked  experimentally  for  self- 
consistency),  »■„  rp,  r'p  (=  ;/>m),  Ep,  F'p(=vp/p), 
and  p'  (~v'p/Ep).  The  primes  denote  an  alternate  method 
of  obtaining  these  parameters.  Comparison  of  rp  with 
rp,  Ep  with  F'p,  and  p with  p'  provided  a critical  check 
of  our  technique.  The  agreement  was  never  worse  than 
15%  and  typically  good  to  within  5%.  Each  datum  point 
plotted  in  Fig.  2 is  an  average  of  the  primed  and  un- 
primed data.  The  sources  of  error  are  uncertainties  in 
changes  of  physical  dimensions  with  pressure  (2%), 

P (5^),  / and  V (3^),  and  p (10^).  It  is  evident  that  rp, 
Fp,  and  r,  all  increase  substantially  with  P up  to  10 
kbar.  In  the  data  shown  in  Fig.  2 the  turnaround  point 
(where  [rp,  F.p\  starts  decreasing  with  P)  is  at  - 10  kbar. 
However,  the  high  R of  the  sample  at  these  /’’s  make 
the  turnover  point  difficult  to  determine  accurately. 

For  comparison  we  also  studied  the  affect  of  P on 
samples  having  intermediate  Fe’s.  VT  either  remained 
constant  or  decreased  slightly  as  P increased  to  10 
kbar.  Again,  no  information  about  the  point  (rp,Ep) 
can  be  extracted  from  intermediate  Fc  samples. 

Vinson  rt  til.  show  that  Monte  Carlo  calculations 
using  /.  as  the  lowest  subsidiary  minimum  provide  good 
agreement  with  the  observations  that  hydrostatic  P 
initially  increases  Ft  and  uniaxial  stress  measurements 
produce  decidedly  different  behavior  for  surfaces  con- 
taining (100)  and  (111)  crystallographic  planes.  Our 


FIG.  2.  Normalized  |>cnk  velocity  (dots)  and  saturated  velocity 
(straight  lines)  versus  normalized  peak  electric  field  as  a 
function  of  P.  The  data  were  extracted  from  I{V)  curves  tor 
three  samples  taken  with  100-nscc  voltage  pulses  at  a prf  of 
50  Hz2.  For  a typical  sample  e^-1.97  xlo7  cm/sec  and  Ep 
-4800  V/cm.  The  curved  line  beginning  at  the  P 0 point  rep- 
resents an  average  of  the  data  appearing  in  Ref.  1.  The  inset 
shows  the  shape  of  a typical  sample,  drawn  to  scale.  The  pre- 
sculptured  dimensions  are  given  in  the  text.  Sculpturing  typi- 
cally reduced  the  cross-sectional  area  by  about  a factor  of  7. 

W.  Czubatyj  and  M.P.  Shaw  206 


hydrostatic  P results  are  in  fair  qualitative  agreement 
with  their  F.p  values  (which  they  obtain,  inappropriately, 
from  Vr)  but  disagree  with  their  results  for  rp,  which 
have  rp  decreasing  smoothly  with  P (Fig.  2).  Further- 
more, our  results  show  that  tlF.JtlP  is  much  larger 
than  they  thought.  Because  of  these  discrepancies,  the 
results  of  their  uniaxial  stress  measurements  must 
also  be  held  as  questionable  since  in  those  cases  the 
samples  were  also  not  shaped  to  eliminate  contact 
effects. 

Without  numerical  calculations,  it  is  difficult  to 
understand  our  data  in  terms  of  the  known  band  struc- 
ture of  GaAs.  The  (vp,Fp)  point  is  observed  to  increase 
very  rapidly  with  P,  faster  than  would  be  expected 
solely  from  an  increase  in  effective  mass  in  the  lowest 
minimum  (D.  It  is  known  that,  with  respect  to  the 
valence-band  maximum,  the  T point  increases  rapidly 
with  P (riF./tlP lr  =-11  xlO'3  eV/kbar)  and  the  X point  de- 
creases slowly  with  P ( rfE/dP\x — 1 xlO"3  eV/ 
kbar). It  is  also  thought  that  the  L point  increases 
with  P,  but  at  a rate  slower  than  the  rate  of  rise  of 
T. 10  If  this  is  the  case,  the  T-L  separation  will  de- 
crease with  P and  hence,  based  on  these  simple  argu- 
ments, (rp,  Fp)  should  decrease.  Clearly  an  appreciation 
of  the  coupling  constants  and  their  P dependence  is  re- 
quired and,  in  view  of  the  recent  findings  of  Aspnes 
cl  rt/.3,  perhaps  the  motion  of  the  conduction-band 
extrema  with  P should  be  reexamined. 

In  conclusion,  two  points  are  evident.  First,  as  pre- 
viously stated,  our  data  were  taken  on  samples  whose 
resistance  increased  exponentially  starting  above  about 
1 kbar.  This  behavior  might  lower  the  mobility  by 
changing  the  amount  of  impurity  scattering  (which  is 
relatively  unimportant  at  300‘K).  Because  of  this  pos- 


sibility, the  changes  in  i pand  v,  with  P might  actually 
be  smaller  than  indicated.  However,  the  measured 
changes  in  F.p  would  be  unaffected  since  they  are  taken 
directly  off  the  raw  data.  For  further  verification,  the 
experiments  should  be  redone  with  sculptured  bulk  or 
shaped  epitaxial  samples  where  /’-induced  carrier 
freeze-out  is  but  a minor  effect.  Second,  the  Monte 
Carlo  calculations  should  be  repeated  in  an  attempt  to 
fit  our  data.  We  have  initiated  such  a program. 

M.P.  Shaw  is  grateful  to  A.  Adams,  W.  Paul,  and 
P.  Vinson  for  useful  conversations.  The  authors  would 
also  like  to  thank  R.  Jeffery  and  D.  Gustafson  for  ex- 
perimental assistance  and  guidance. 


•Research  supported  by  the  Office  of  Naval  Research  uiylcr 
contract  No.  N-0014-75-C-0399.  , 

1 P. J . Vinson,  C.  Pickering,  A.R.  Adams,  W.  Fawcett,  and 
G.D.  Pitt,  Proc.  tilth  Int.  Conf.  Phys.  Semiconductors, 
Rome,  1976  (unpublished).  ‘ 

2Sec,  e.g.  , P.R.  Solomon,  M.P.  Shaw,  ILL.  Grubin,  and 
R.  Kaul,  IF1K  Trans.  Klectron  Devices  ED-22,  (1975). 

’D.E.  Aspnes,  G.G.  Olson,  and  D.W.  I.ynch,  Phys.  Rev. 

Lett.  37,  766  (1976). 

*C.  Pickering,  A.R.  Adams,  G.l).  Pitt,  andM.K.R.  Vyas, 

J.  Phys.  C 8,  129  (1975). 

’A.R.  Hutson,  A.  Jayaraman,  A.G.  Chynowcth,  A.S.  Coricll, 
and  W.C.  Feldman,  Phys.  Rev.  Lett.  14,  639  (1965). 

®M.  P.  Shaw,  P.R.  Solomon,  and  ILL.  Grubin,  IBM  J.  Res. 
I)cv.  13,  587  '1969). 

*11.1..  Grubin,  M.P.  Shaw,  and  P.R.  Solomon,  IEEE  Trans. 
Electron  Devices  ED-20,  63  (1973). 

8P.  Guctin  and  G.  Schroder,  Phys.  Rev.  B 5,  3979  (1972). 

G . I) . Pitt  and  J.  I.ees,  Proc.  of  Lcs  Proprieties  Physiques 
tics  Snliilcs  Sons  Pression  (Centre  National  de  la  Recherche 
Scientiiique.  Paris,  1970),  p.  225. 

10D.  L.  Camphausen,  G.A.N.  Connell,  and  W.  Paul,  Phys. 

Rev.  Lett.  26,  184  (1976). 


207 


Applied  Phyiics  Letter!,  Vol.  30,  No.  4.  15  February  1977 


207 


